General form of quantum evolution 
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We propose a complete treatment of a local in time dynamics of open quantum systems. In 
this approach Markovian evolution turns out to be a special case of a general non-Markovian one. 
We provide a general representation of the local generator which generalizes well known Lindblad 
representation for the Markovian dynamics. It shows that the structure of non-Markovian generators 
is highly intricate and the problem of their classification is still open. Simple examples illustrate 
our approach. 
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Any realistic quantum system inevitably interacts with 
its environment, therefore, the theory of open quantum 
systems and their dynamical features is of particular im- 
portance p], Q. Actually, this problem attracts nowa- 
days increasing attention due to the growing interest in 
controlling quantum systems and applications in modern 
quantum technologies such as quantum communication, 
cryptography and computation Q. 

For several decades the popular Markovian approxima- 
tion which does not take into account memory effects was 
successfully studied and applied in a variety of problems 
P, 0. However, recent investigations in quantum infor- 
mation and recent technological progress calls for truly 
non-Markovian approach. In the last few years many 
analytical methods and numerical techniques have been 
developed to treat non-Markovian processes in quantum 
optics, solid state physics and quantum information 0- 
Il7| . Moreover, several measures of non-Markovianity 
were proposed and intensively studied [l8l - [2^ |. 

The most general form of local in time Master Equa- 
tion reads as follows 



dp(t) 
dt 



C(t,t )p(t) , p(t Q ) = p 



(1) 



where £(i, to) is a local generator which depends not only 
upon the current time 't' but in principle it might depend 
upon the initial point 't '. It is clear that dependence on 
'to' introduces an effective memory. The system does 
remember when the evolution begun. We call the evo- 
lution governed by ((T|) Markovian if and only if £(t,t ) 
does not depend on 'to'- Otherwise the evolution is non- 
Markovian. It is clear that £(t,to) is defined for t > to 
only. Note however, that in the Markovian case £m(£) 
is defined in principle for all t G (— oo, oo). Any solution 
to ((T|) gives rise to the dynamical map A(t,t ) defined 
by p(t) = A(t,to)po- Clearly A(t,to) itself satisfies the 
following equation 



dt 



A(t, t ) = C(t, t )A(i, t ) , A(to, t ) = 1, (2) 



where II denotes an identity map. A particular class of 
£(t,to) is provided by the homogeneous generators, i.e. 
when £ depends on 't — to' only. Hence, the evolution 



is governed by the 1-paramctcr family £(r) defined for 
t > 0. It is clear that in this case the evolution is ho- 
mogeneous as well, that is A(t + T,t + T) = A(t, to) 
for arbitrary T, and hence one may define a 1-paramcter 
family of dynamical maps A(t) := A(t, 0). Actually, one 
usually fixes to = from the very beginning and consid- 
ers 



j t A(t) = £(t)A(t) , A(0) = 11 



(3) 



We stress that it could be done only in the homogeneous 
case and usually it is referred as a time convolutionlcss 
(TCL) approach ^MM- 

A solution A(t,t ) to @ is defined by the following 
formula 

A(t, t ) = T exp £(r, t )d^j , (4) 

where T stands for the chronological operator. In the 
homogeneous case it simplifies to 

A(Mo) = Texpf / C(T)dr) , (5) 



which clearly shows that A(t, to) depends on H — to'- We 
stress that the formula (|4|) has only a formal character 
since in general the evaluation of T-product is not feasi- 
ble. Recall, that this formula simplifies if C(t, to) defines 
mutually commuting family, i.e. [C(t, to), C(u, to)] = 
for all t, u > to- In this case T-product drops out from 
©• 

A solution A(t, to) defines a legitimate quantum dy- 
namics if and only if A(t, to) is completely positive and 
trace preserving (CPT) for all t > to- Now comes the 
natural question: how to characterize the properties of 
£(t,to) which guarantee that A(f, to) corresponds to the 
legitimate quantum dynamics. These conditions are well 
known in the Markovian case: a solution to 



dt 



A(t, t ) = £ M (t)A(t, t ) , A(t , t ) = 1 , (6) 



is CPT if and only if the time dependent generator has 
the following Lindblad representation [26l - l28{ 

C MP = -i[H, p] + 5>„ (VcpVZ - \{V} t V a ,p}) , (7) 
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whore H = H(t) stands for the effective time-dependent 
Hamiltonian and V a = V a (t) are time-dependent Lind- 
blad (or noise) operators. The time dependent coeffi- 
cients 7 Q satisfy 7 Q (t) > and encode the information 
about dissipation and/or decoherence of the system. 

Let us observe that a family A(t, to) of CPT maps may 
1)0 represented by 

A(t,t ) = e z ^ , (8) 

where Z[t, to) has a Lindblad representation for all t > to- 
The price we pay for this simple representation is that 
Z(t,to) might be highly singular. It is clear that formally 
Z(t,to) is defined as a logarithm of A(t, to) and hence it 
is not uniquely defined (log has an infinite number of 
branches). Moreover, one always meets problems when 
A(t, to) possesses eigenvalues belonging to the cut of log, 
cf. discussion in [18|. For example a CPT map Ap = 
o z pa z cannot be represented by A = e z . Note, however, 
that A may be considered as a limit of A(t) = e z ^\ 
with Z(t) = — log(cos t)Lo, and Lop = <J z pVz — p is a 
legitimate Lindblad generator for t £_[0,7r/2). One has 
A = lim t ^^-/2 A(t) (see discussion in jl7j|). 

Note, that condition A(io,io) = 1 is equivalent to 
Z{to, to) = which is guarantied by 

Z(t,t Q ) = f X{u,t )du , (9) 

and hence the solution has the following form 

A(t,t ) = cxp Qf X(T,to)dT^j . (10) 

Note, that contrary to (fj| the above formula does not 
contain chronological T-product. The corresponding gen- 
erator C(t,t ) is defined by [ItJ 

£(t,t Q ) = ^A(Mo) • A(Mq)- 1 , (11) 
at 

where A(t,to) _1 = e^ z< - t ' t " s> denotes the inverse of 
A(t,t ). Note, that A(t,t )~ 1 is not completely positive, 
hence can not describe quantum evolution backwards in 
time, unless A(t, to) is unitary or anti-unitary. Now, to 
compute dA(t,to)/dt one uses well known formula [29j 

^e A ^ = f e sA ^A{t)e^ A ^ds , (12) 
dt Jo 

where A(t) is an arbitrary (difFerentiable) family of op- 
erators, and A = dA/dt. Hence 

* zfcto) =£f t ,to)e z ^ , (13) 
dt 

where 

£(t,to)= [ e sZ ^ to) X(t,t )e- sZ( - t ' to ' ) ds . (14) 
Jo 



This is the main result of our Letter. It proves that each 
legitimate generator £(t,to) of quantum evolution gov- 
erned by the Master Equation (flj has the form defined 
by (fT4"]) . where Z(t,to) has a Lindblad representation for 
each t > to, and X(t, to) is defined in (|9]). Hence the con- 
struction of a legitimate generator is pretty simple: each 
family of Lindblad operators Z(t,to), with Z(to,to) = 0, 
gives rise via (|14[) to the corresponding prescription for 
£(t,to)- Nevertheless, the formula (fl4|) is highly nontriv- 
ial and the computation of £(t,to) out of Z(t,to) might 
be highly complicated. This is the price we pay for the 
simple representation of evolution (|10p . Hence, we have a 
kind of complementarity: either one uses T-product for- 
mula (|4]) with relatively simple generator or one avoids 
T-product in ([10j) but uses highly nontrivial generator 
(TH| . The advantage of our approach is that one knows 
how to construct generator (in practice it might be com- 
plicated) giving rise to the legitimate quantum dynamics. 

Let us observe that in the special case when X(t, to) 
mutually commute, i.e. [X(t,to),X(u,to)] = for all 
t, u > to, the formula (|14l) reduces to C{t, to) = X(t, to) . 
Hence, a commuting family C(t,to) defines a legitimate 
generator if and only if Z(t,to) = J t C(u,to) du has a 
Lindblad representation for all t > to- In the noncom- 
mutative case this simple criterion is no longer true. 

The characteristic feature of the Markovian evolution 
governed by ([6]) is that A(t, to) satisfies local composition 
law 

A(t,s)-A(s,to) =A(Mo) , (15) 

for t > s > t . Actually, this property is guarantied by 
the intricate action of T-product in the formula (U]). Now, 
changing the representation from (|4|) into (| 10[) the valid- 
ity of composition law is no longer visible. The formula 
(fl~5|) implies 

e Z(t,u) , e Z(u,t ) _ e Z(t,ta) ^ 

for t > s > to Clearly, in the commutative one has simply 

Z(t,u) + Z(u,t ) = Z(t,t Q ) . (17) 

Note however that when Z(t,to) do not commute, the 
Baker-Campbcll-Hausdorff formula e A e B = e c , with 

C = A+B+ l -[A,B] + ^([A, [A, B]] - [B, [A,B}})+..., 

provides highly nontrivial condition upon Z(t,to). It 
shows that knowing legitimate Z(t, to) one does not 
know immediately wether the corresponding dynamics is 
Markovian or not. Only applying ([T4"]) one can check 
wether £(t,to) does, or does not, depend on 'to' and 
hence infer about Markovianity. This problem simpli- 
fies in the homogeneous case: now the evolution is never 
Markovian unless Z(t,to) = (t — to)Lo, i.e. £(t,to) = Lq, 
where Lo denotes the time independent Markovian gen- 
erator. 



3 



It is clear that for a general family of Lindblad gen- 
erators Z(t,t ) computation of £(t, to) via (fLE)) is rather 
untractable. However, us usual, there is a class of Z(t, to) 
for which the general problem simplifies considerably. 
Consider a special class of Z(t, to) defined via © by the 
following family 

X(t,to) = a 1 (t,t )Ci + ... + a N (t,to)C N , (18) 

where Ci, . . . , £n are time independent Lindblad gener- 
ators. One has 

Z(t, t ) = Ai(t, t )A + . . . + A N (t, t )C N , (19) 

where Ak(t, to) = f* g ak(u,to)du. Now, Z(t,to) has a 
Lindblad representation iff Ak(t,to) > 0. Let us observe 
that if C\, . . . , Cn close a Lie algebra, i.e. [£j,£j] = 
Sfe=l c ij£k, then using well known Lie algebraic methods 
one can easily compute £(t,to) out of ([H]) and gets 

C(t, to) = b^t, t )£i + ... + b N (t, t )C N , (20) 

where the functions bk{t,to) are uniquely defined by 
Qk{t,to) and the structure constants c*j. Actually, any 
set {£i, . . . ,Cn} of Lindblad generators may be always 
completed to close a Lie algebra. It follows from the 
fact that a set of Lindblad generators belong to the 
Lie algebra corresponding to the Lie group of linear 
maps preserving hcrmiticity. Note, that if [Ci,Cj\ = 0, 
i.e. the corresponding Lie algebra is commutative, then 
bk(t,h) = a k (t,t ), that is, £{t,t ) = X(t,t ). 

Example 1 (Commutative case) Consider the fol- 
lowing pure decoherence model defined by the follow- 
ing time dependent Hamiltonian H(t) = Hn(t) + 
Hs(t) + Hsn(t), where H R {t) is the reservoir Hamil- 
tonian, H s (t) = ^„e„(t)?„ (P„ = \n)(n\) the system 
Hamiltonian and 

H SR (t) = Y^Pn®B n (t) (21) 

n 

the interaction part, B n — B^ being reservoirs operators. 
The initial product state p®ui evolves according to the 
unitary evolution U(t,to)(p®uj)U(t,to)' < and by partial 
tracing with respect to the reservoir degrees of freedom 
one finds for the evolved system density matrix 

p(t) = A(t, to)p = ^ C mn(t, to)P m pPn , (22) 

n , m 

where c mn (t,t ) = Tr[U m (t, t ) u U n (t, t ) t ], with 
U n (t,t ) = TexpH Jl Yn(r)dr], and Y n (r) = e n (r)I R + 
Hr(t)+B„(t) being time dependent reservoir operators. 
Note that the matrix c mn (t,to) is semi-positive definite 
and hence (|22|) defines the Kraus representation of the 
completely positive map A(t,t ). Note that A(t,t ) de- 
fines a commutative family of maps and hence one easily 
finds for the corresponding generator 

£(i, t ) p = ^ a mn{t, t )P m pPn , (23) 



where the functions a mn (t,to) are defined by 
a mn (t,t ) = c mn (t,to)/c mn (t,to). Note that if 
dimHs = 2, then c\\ = C22 = 1 , and C12 = 7 
with |7| < 1. One easily finds for the local generator 

C(t,t )p = ibi(t,t )[a z ,p] - b 2 (t,to)[cr z pa z - p] , (24) 

with bi = Im(7/27) and b 2 = Re(7/27). Note that this 
dynamics is homogeneous if and only if the Hamiltonian 
of S + R is time independent. 

Example 2 (Noncommutative case) Let us consider 
a simple example of exactly solvable dynamics of 2-level 
system defined by the following homogenous family of 
operators X(t) := X(t,0) 

X(t) = a 1 (t)£ 1 +a 2 (t)£ 2 , t>0, (25) 

where the Markovian generators C\ , £2 arc defined by 

dp = a + p(j~ - i{cr~o-+,p} , 

C 2 p = (J-pa + -i{cr+(T-,p} , 

and <j + = |1)(2| , a~ = |2)(1| are the standard raising 
and lowering qubit operators ({|1), |2)} denotes an or- 
thonormal basis in the qubit Hilbert space). Since C\ 
and £2 do not commute X(t) defines a noncommutative 
family. Clearly, one may add to X(t) a commutative 
part (IM1) which does commute with X(t) and hence do 
not change qualitative features of dynamic. For simplic- 
ity we consider only simplified version which is essential 
for non-commutativity. 

The time dependent parameters a\{t) and a 2 (t) are 
arbitrary but real. Following our construction one has 

Z(t) = A 1 (t)C 1 +A 2 (t)C 2 , (26) 

where Ak{t) = L ak(u)du. Hence, the formula A(t) = 
e z W defines CPT map for all f > if and only if 

A x {t) > , A 2 {t) > . (27) 

Now, let us apply our basic formula (|14p to find the cor- 
responding generator £(t). Observing that £\ and £ 2 
close a Lie algebra [£i, £ 2 \ = £\ — £ 2 , one easily finds 

£(t)=b 1 (t)£ 1 +b 2 (t)£ 2 , (28) 

where 

M*) = oiW + /(*) > b 2 (t)=a 2 (t)-f(t) , (29) 
and the time dependent function f(t) reads as follows 

where the Wronskian W(t) = Ax{t)a 2 {t) - A 2 {t)a x {t), 
and A(t) = A x {t) + A 2 (t). Note, that f(t) = (for all 



t > 0) if and only if W(t) = 0, i.e. functions ai(t) and 
a 2 (t) are linearly dependent. If this is the case one has 
a 2 (t) = Aai(i), and X{t) = ai(t){Ci + \£ 2 ) defines a 
commutative family. In this case one has bk{t) = a,k{t) 
and hence £{t) = X{t). In the general noncommutative 
case one has for the integral 

/ C(r)dT = Bi(t)& + B 2 {t)C 2 , (31) 
Jo 

with Btit) = A 1 (t) + F(t) , B 2 (t) = A 2 (t) - F(t) , and 
F(t) = f*f(u)du. Note, that Bi(t) + B 2 (t) = A(t) > 
0. However, contrary to the commutative case, there 
is no need that both Bi{t) and B 2 (t) are positive. It 
shows that integral J £(r)dr needs not have a Lindblad 
representation. Note, that 



£(t) = X(t) + f(t)[£ 1 -£ 2 ] , 



(32) 



which clearly shows that the last term '/(i)[£i — £ 2 \ 
destroys the Lindblad structure of J* £{t)cLt. It proves 
the intricate action of T-product: 



Texp (j [X(r) + /(r)[£ a - £ 2 }}d^j 



(33) 



that is, chronological product simply washes out the un- 
wanted term l f(t)[£i — £ 2 }\ Eventually, one easily shows 
(using standard algebraic methods, e.g. [3(|) that 



lni/i(t)£i _ lnu 2 (t)C 2 



(34) 



where (skipping time dependence) 



.4 



-A 



A 2 



l> 2 



Ax + A 2 e J 



(35) 



Note that v\v 2 — e A . One has Vk(t) > 1, and hence 
P4p gives another representation of dynamical map as a 
composition of two completely positive maps generated 
by £\ and £ 2 . 



In conclusion, we proposed a complete treatment of a 
local in time dynamics of open quantum systems based 
on the Master Equation (QJ . We provided a general rep- 
resentation of the local generator - formula (|14j) - which 
generalizes well known Lindblad representation for the 
Markovian dynamics. We stress that any local generator 
£(t,to) may be constructed via (|T4j) by a suitable choice 
of the Lindblad family Z(t,to). However, the problem of 
necessary and sufficient condition for £(t,to) which guar- 
antee that A(t,to) is CPT is still open. Only, if £(t,to) 
defines a commutative family, these conditions reduce to 
a simple requirement that J t £(u,to)du has a Lindblad 
form for t > to. 
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